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We show that, in the context of brane-world scenarios with low tension r = / 4 , massive brane 
fluctuations are natural dark matter candidates. We calculate the present abundances for both 
hot(warm) and cold branons in terms of the branon mass M and the tension scale /. The results 
are compared with the current experimental bounds on these parameters. We also study the 
prospects for their detection in direct search experiments and comment on their characteristic 
signals in the indirect ones. 
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PACS numbers: 95.35.+d, ll.25.-w, ll.10.Kk 

One of the most important open problems in astro- 
physics and cosmology is to identify the nature of dark 
matter. It has been known for a long time that the lu- 
minous matter observed in spiral galaxies is insufficient 
to explain their rotation curves. The existence of dark 
halos was proposed as a possible solution for the discrep- 
ancy (see [l( and references therein), although at present, 
numerical simulations of the formation of such halos ap- 
pear inconsistent with observations. On the other hand, 
different estimations of the total matter density in the 
universe from large scale motions, virial masses or clus- 
ter abundances, and the more recent Type la super- 
novae and CMB anisotropics observations agree in a value 
tt M = 0.27 ± 0.04 which is much larger than the 

value of the total luminous mass density in the universe 
Qiumh — 0.006 — 0.002. In addition, the big bang nucle- 
osynthesis (BBN) and WMAP results for the total bary- 
onic content n B h 2 = 0.0224 ± 0.0009 imply that most of 
the matter in the universe is dark and non-baryonic (see 
HIS and references therein) . 

A possible explanation of this puzzle is that the dom- 
inant component of dark matter consists of some non- 
relativistic (cold) stable and weakly interacting massive 
particles (WIMP) which decoupled from radiation early 
enough so that their relic abundances are important to- 
day. The possibility that the universe is dominated by 
hot dark matter seems to conflict with numerical simu- 
lations of structure formation. Thus, the only potential 
candidates within the known particles would be massive 
neutrinos. However a detailed analysis has excluded the 
three light and even one additional heavy fourth gener- 
ation of Majorana or Dirac neutrinos. This fact has led 
to the search of cold dark matter candidates beyond the 
Standard Model (SM) 0. 

There are two main such particles studied in the litera- 
ture. On one hand axions which strictly speaking cannot 
be considered as WIMP's since they are very light and 
produced non-thermally. On the other hand we have the 
lightest supersymmetric particle, which can be identified 
with a neutralino in most of the supersymmetric models. 
The latter is probably the most studied and best theoret- 



ically motivated dark matter candidate [4| . However the 
large number of free parameters in supersymmetric the- 
ories make their predictions extremely model dependent. 
More recently, the existence of large extra dimensions has 
been proposed as a new setting for a possible solution 
to the hierarchy problem In this scenario, the SM 
fields are forced to live on a three-dimensional hypersur- 
face (brane) whereas gravity is able to propagate on the 
higher D = 4 + N dimensional bulk space. In this Brane 
World scenario (BWS) the fundamental scale of gravity 
is not the Planck scale Mp but another scale Mr>. Al- 
though in the original ADD model Q, Mr> is supposed 
to be not too much larger than the electroweak scale, re- 
cently brane cosmology models have arised in which Mrj 
is much larger than the TeV (||| and references therein). 
In this work we will consider the case of a general BWS 
without assuming any particular value for Mr). 

In these models gravitons propagating through the 
bulk space give rise to a Kaluza-Klcin (KK) tower of mas- 
sive gravitons on the brane. These KK gravitons couple 
to the energy- momentum tensor of the SM fields Tg V M and 
could be produced under the appropriate circumstances 
as real or virtual particles. Another important effect that 
is expected in the BWS is the presence of brane fluctu- 
ations since rigid objects do not exist in relativistic the- 
ories. In other words the brane should have some finite 
tension t = / 4 . When these oscillations are taken into 
account two new effects appear 0. First of all, we have 
to introduce new fields, which for a homogeneous extra 
space, essentially represent the position of the brane in 
the bulk space (x M , y a ~ n a (x)/ f 2 ). The n a (x) fields are 
the Goldstone bosons (GB) corresponding to the spon- 
taneous symmetry breaking (SSB) of the translational 
invariance produced by the presence of the brane (bra- 
nons) . It has been shown Q that when these branons are 
properly taken into account, the coupling of the SM par- 
ticles to any bulk field is exponentially suppressed by a 
factor exp(-M|- K M| ) /(87r 2 / 4 )), where M KK is the mass 
of the corresponding KK mode. As a consequence, if the 
tension scale / is much smaller than the fundamental 
scale Mr), i.e. / <C Md, the KK modes decouple from 
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the SM particles. Therefore, for flexible enough branes, 
the only relevant degrees of freedom at low energies in the 
BWS are the SM particles and the branons. Similarly to 
other GB's, branons are expected to be nearly massless 
and weakly interacting at low energies. Nevertheless, in 
general, translational invariancc in the extra dimensions 
is not necessarily an exact symmetry and some branon 
mass M is expected from such explicit symmetry break- 
ing as shown in 0, E| . This is similar to what happens 
to pions which are the GB corresponding to the SSB of 
the chiral symmetry of low-energy strong interactions. 
As gravitons do, branons couple to Ttt M , however in this 
case, the lowest order effective Lagrangian is 0: 

+ ^(W^d^-MWg^T^ (1) 

We see that branons always interact by pairs, they are 
stable and difficult to detect, since their interactions are 
suppressed by the tension scale /, and they are expected 
to be massive. Thus we arrive to the conclusion that the 
massive oscillations of the brane are natural candidates 
to dark matter in the BWS where / <C Mb. The dark 
matter problem has been considered in this scenario from 
a different point of view in llOl and also in models with 
universal extra dimensions 

In order to calculate the thermal relic branon abun- 
dance, we will use the standard techniques given in 
in two limiting cases, either relativistic (hot) or non- 
relativistic (cold) branons at decoupling. The evolution 
of the number density n a of branons 7r a , a = 1 , . . . , TV 
interacting with SM particles in an expanding universe 
is given by the Boltzmann equation: 



dn a 
~dt 
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where a a = J2x °'{' Ka K a —> X) is the total annihilation 
cross section of branons into SM particles X summed over 
all final states. The —3Hn a term, with H the Hubble 
parameter, takes into account the dilution of the num- 
ber density due to the universe expansion. These are the 
only terms which could change their number density to 
the leading order. In fact, branons do not decay into 
other particles and since they interact always by pairs 
the conversions like ir a X — > ir a Y do not change their 
number. Notice that we are considering the low-energy 
effective Lagrangian in Q and assuming for simplicity 
that all the branons are degenerate. Accordingly, each 
branon species evolves independently, and in the follow- 
ing we will drop the a index. The total branon density 
will be just N times that of a single branon. The (<tav) 
term denotes the thermal average of the total annihila- 
tion cross section times the relative velocity. From Q , it 
includes, to leading order, annihilations into all the SM 
particle-antiparticle pairs. If the universe temperature is 



above the QCD phase transition (T > T c ), we consider 
annihilations into quark-antiquark and gluons pairs. If 
T < T c we include annihilations into light hadrons. For 
the sake of definiteness we have taken a critical temper- 
ature T c ~ 170 MeV and a Higgs mass tuh — 125 GeV, 
although the final results are not very sensitive to the 
concrete value of these parameters. 

Defining the new variable x — M/T, one finds from 
(0) that below the freeze-out temperature Xf for which 
the annihilation rate Ta = n eq {cr A v) equals the ex- 
pansion rate H, branons are decoupled from the ther- 
mal bath and their abundance remains frozen relative 
to the entropy density. We will denote g e ff(T) and 
h e ff(T) the effective relativistic degrees of freedom for 
energy and entropy densities respectively at tempera- 
ture T. In order to calculate the decoupling temper- 
ature for hot relics, it is a good approximation to use 
the condition YA{xf) — H(xf). From the Friedmann 
equation in a radiation dominated universe we have: 
H(x) ~ 1.67 gl f / f M 2 /{x 2 M P ) On the other hand, 
expanding Ta(x) for x <C 3 and neglecting M, we 
find: T\{T) ~ 16tt 9 T 9 /(297675C(3)/ 8 ) for photons and 
T A (T) — Y A (T)/A for neutrinos. For massive particles 
we cannot give closed expressions. Once we know x/, 
the corresponding fraction of energy density today in the 
form of relics is given by: 
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h eff(xf) eV 
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For cold relics the calculation of the decoupling temper- 
ature is more complicated. The well-known result is: 



x f 



In 



0.038 c(c+2)M P M(a A v) 
T/2 372 
9 eff X f 
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where c ~ 0.5 is obtained from the numerical resolution 
of the Boltzmann equation. The above equation can be 
solved iteratively. The matter density can be written as: 



f2 Brh 
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where we have expanded (<jav) i n powers of x^ 1 as 
(cr A v) = J2^=o c n x ~ n ■ I n the case of photons, the first 
non-vanishing coefficient is — 68A/ 6 /(15/ 8 7r 2 ) and 



for massless neutrinos 



- P 



(d— wave annihilation) 



whereas for non-conformal matter we also have s— and 
p— wave annihilation, (co,ci 7^ 0). The corresponding 
expressions are more complicated and will be given else- 
where. We have performed all the expansions up to 
0(x~ 2 ). Coannihilation effects are absent in this case 
since there are no slightly heavier particles which even- 
tually could decay into the lightest branon. Also, in order 
to avoid the problems of the Taylor expansion near SM 
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thresholds, we have taken branon masses sufficiently sep- 
arated from SM particles masses where the usual treat- 
ment is adequate |l2| . Such treatment is known to intro- 
duce errors of the order of 10% in the relic abundances. 



10 



> 10' 

a 



10 -, 



■111 


■I 




I I i I 



10 



10 



-5 



-1 



10 

M (GeV) 



10' 



10 



FIG. 1: Relic abundance for a model with one branon N — 1. 
The line on the left is the Q Br h 2 = 0.129-0.095 curve for hot- 
warm relics. The line on the right corresponds to cold relics. 
The lower striped area is the estimated excluded region by 
single-photon processes at LEP-II [9j and the upper area is 
also excluded by branon overproduction. 

Branons could be responsible for the observed cosmo- 
logical dark matter density provided Q Br h 2 — 0.129 — 
0.095 at the 95% C.L. which corresponds to fls r = 

0. 23 ± 0.08 and h = 0.79 - 0.65 0. In FIG. 1, we 
have plotted these curves for hot and cold branons in 
the / — M plane for one single branon. For N types 
of branons the corresponding abundances are simply N 
times larger. In fact, the contribution of branons to g e ff 
in J5J is negligible in the case of cold relics. For hot 
relics such contribution in © has been taken into ac- 
count, although it is very small in the interesting re- 
gions. Concerning the freeze-out temperature, the re- 
sults on the cold dark matter curve range from xt ~ 9 
for M = 10~ 5 GeV to x f ~ 31 for M = 10 6 GeV. In 
the hot relics case, a very good approximation is given 
by log(7>/GeV) ~ (8/7) log(//GeV) - 3.2. 

Pure hot dark matter models are disfavored at present 
because relativistic matter free streams from overdense 
into underdense regions preventing structures from grow- 
ing below the so called free-streaming scale given by 0] : 
X FS ~ 0.2(fi Bl .^ 2 ) 1 / 3 (keV/M) 4 / 3 Mpc. In the hot bra- 
nons curve in FIG. 1, the masses in the allowed region 
are in the range M = 85 — 177 eV, which corresponds to 
Xfs — 2.4 — 1.0 Mpc. Such scales are much smaller than 
those in neutrino dark matter models and, in addition, 
since these branons decouple much earlier than neutri- 
nos do, their corresponding temperatures are also lower, 

1. e. they could be considered rather as warm dark matter 



(WDM) candidates from the point of view of struc- 
ture formation. However, a WDM dominated universe 
seems to be also disfavored by the recent observations 
of WMAP 0. On the other hand, the presence of rela- 
tivistic branons during BBN could change the expansion 
rate of the universe, spoiling the predictions of the light 
elements abundances [l2j. However, the BBN limits on 
the number of branons contributing to g e f / are not very 
constraining. In fact, imposin g th e conservative bound 
9eff{T B BN ~ lMeV) < 12.5 0, we find N < 18 for 
branons which decoupled before the QCD phase transi- 
tion, corresponding to / > 60 GeV. 

Let us clarify the main assumptions used so far in the 
work. We are considering models with / -C Mjj , and also 
assuming that the evolution of the universe is standard 
up to a temperature around /. Indeed, this is the case 
of realistic brane cosmology models |6J. Therefore, we 
have taken the conservative bound / > 1 MeV, so that 
BBN is not affected. Moreover the effective Lagrangian 
JU is only valid at low energies relative to /. We have 
checked that our calculations are consistent with these as- 
sumptions since the decoupling temperatures are always 
smaller than / in the allowed regions in FIG. 1. 

Brane fluctuations could be, not only candidates for 
the cosmological dark matter, but also they could make 
up the galactic halo and explain the local dynamics. In 
such case, they could be detected in direct search exper- 
iments from the energy transfer in elastic collisions with 
nuclei of a suitable target. The appropriate quantity to 
be compared with the experimental results is not the elas- 
tic branon- nucleus cross section cr, but the differential 
cross section per nucleon at zero-momentum transfer a n , 
which is defined by 0: 



da _a n A 2 F 2 (\q\ 



d\q\ 



2, .2 



(6) 



where [i = Mm/ '(M + m), F(\q\) is a nuclear form factor 
with the normalization F(0) = 1, m ~ 939 MeV is the 
nucleon mass, v is the relative velocity and A is the mass 
number of the nucleus. In the limit in which the momen- 
tum transfer goes to zero, we can consider the nucleons 
as pointlike particles. In this case, it is possible to cal- 
culate the branon-nucleon cross section er„ from 0) just 
considering the nucleon as a Dirac fcrmion of mass m: 



9M 2 toV 
64tt/ 8 



(7) 



In fact, this quantity does not depend on the type of par- 
ticle which couples to the branon, but only on its mass. 
This can be seen from Q since in this limit, branons only 
couple to the T 00 component. 

The results of our analysis are shown in FIG. 2. Lines 
of constant / with 50 GeV separation are shown for ref- 
erence. The area on the left of the flB r h 2 =0.129 — 0.095 
curves is excluded by branon overproduction, but the 
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FIG. 2: Elastic branon-nucleon cross section a n in terms 
of the branon mass. The two thick lines correspond to the 
n Br h 2 = 0.129 - 0.095 curve for cold branons in Fig.l with 
N — 1 (left) and N = 7 (right). The shaded areas on the 
left are the previous LEP-II exclusion regions 0], also for 
N = 1,7. The solid lines correspond to the current limits 
on spin-independent cross section from direct detection ex- 
periments: ZEPLIN1 0, DAMA [3 and EDELWEISS [T^J. 
The discontinuous lines are the projected limits for: CRESST 
Tvji, GEDEON |3, CDMS J3, ZEPLIN4 |H and GENIUS 
|21j . (limits obtained from (22j|1. 

right portion is compatible with observations and will be 
explored in future experiments. Such region corresponds 
to / > 120 GeV and M > 40 GeV. 

Another interesting possibility is to detect branons in- 
directly. Their annihilations in the galactic halo can 
give rise to pairs of photons or e + e~ which could be de- 
tected by 7— rays telescopes such as MAGIC or GLAST 
or antimatter detectors (AMS). Annihilation of branons 
trapped in the center of the sun or the earth can give 
rise to high-energy neutrinos which could be detectable 
by high-energy neutrino telescopes such as AMANDA, 
IceCube or ANTARES (see for example 11]). Because 
annihilations of non-relativistic branons into conformal 
matter are d— wave suppressed, the most relevant contri- 
bution will come from the secondary leptonic decays of 
ultra-relativistic Z or W bosons (in the case M 3> Mz). 
These processes will be characterized by the presence of 
peaks around one half of the branon mass in the leptonic 
or neutrino spectra. In the case of photons, softer peaks 
will be present at lower energies and therefore their de- 
tection will be more difficult. The hadronic decays will 
give rise to relatively smeared spectra at lower energies. 
Detailed results will be presented elsewhere. 

Throughout the paper we have assumed that branons 



were in thermal equilibrium with radiation at some point 
in the history of the universe. If this is not the case, 
branons could still be produced non-thermally, very much 
in the same way as axionic dark matter |2.'"j| . In fact for 
very light branons, the energy density produced by this 
mechanism could be cosmologically important. 

In conclusion, we have proposed branons as natural 
dark matter candidates in the BWS with low tension. 
Our results show that in a certain range of the parame- 
ters / and M, their relic abundances could explain the 
missing mass problem, and that such parameters region 
will be explored in future direct detection experiments. 
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